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Abstract. We analyze the Majorana equation in the limit where the particle is at 
rest. We show that several counterintuitive features, absent in the rest limit of the Dirac 
equation, do appear. Among them, Dirac-like positive energy solutions that turn into 
negative energy ones by free evolution, or nonstandard oscillations and interference 
between real and imaginary spinor components for complex solutions. We also study 
the ultrarclativistic limit, showing that the Majorana and Dirac equations mutually 
converge. Furthermore, we propose a physical implementation in trapped ions. 
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1. Introduction 

Quantum simulators are quantum controllable devices aiming at the reproduction of 
complex quantum system dynamics Among their promised advantages are an 
exponential speedup in the computational time with respect to classical computers 
for simulating quantum systems, and the possibility to analyze systems difficult 
to experiment with in the lab [2j [3]. Many physical systems have already been 
emulated in a quantum simulator. Among them, spin models [U [5], quantum phase 
transitions [6J, quantum chemistry [7J, spin statistics [8], and relativistic systems 
[SI QUI EED EE21 E31 Ell E51 EES] • Recently, a theoretical proposal for the quantum simulation 
of the Majorana equation and unphysical operations in trapped ions was proposed [T7] . 

The Majorana equation is a non-Hamiltonian relativistic wave equation that takes 
the form [18] 

ihtftip = mcipc, (1) 

where $ = j^d^ and 7 M are the Dirac matrices [15] , and ip c is the charge conjugate of 
ip. More explicitly, for 1+1 dimensions, the Majorana equation is written as 

ih^doip + ih^ l d x ip = mcip c , (2) 

where the 2x2 matrices 

and 

i = 7 1 7V- (4) 
Notice that, as opposed to the Dirac equation, this is a non-Hamiltonian equation due 
to the simultaneous presence of ip and ip c . Thus, this equation contains counterintuitive 
features, like the nonexistence of energy eigenstates. In Ref. [17], it was shown that 
the Majorana equation can be mapped to a higher-dimensional Dirac equation that has 
a Hamiltonian evolution, such that the Majorana dynamics may be implemented in a 
controllable quantum optics setup as trapped ions [T7j . 

In this work, we go a step further and analyze the strict nonrelativistic limit in which 
the Majorana particle is at rest. Thus, in Section [2] we consider the limit p/mc <C 1, 
where all moments of the particle momentum p are negligible. We will show that, 
surprisingly, novel features emerge, even in the rest limit of the particle p ~ 0: if we 
begin with a Dirac-like spinor with positive energy, the Majorana equation will mix its 
dynamics with the equivalent to the Dirac equation negative energy branch, and it will 
produce nonstandard interference between real and imaginary components. Thus, the 
rest limit of the Majorana equation behaves in a totally different way to the equivalent 
limit of the Dirac equation. Indeed, this is the limit where both equations are more 
different. In Section [3j we show how to implement this dynamics in a trapped-ion system 
and, in Section |4j we present our conclusions. Finally, for the sake of completeness, we 
include in the Appendix the exact solution of the Dirac and Majorana equations and 
its mutual convergence in the ultrarelativistic limit in which mc <C p. 
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2. The rest limit of the 1+1 Majorana equation 

The dynamics of a general spinor for the 1+1-dimension Majorana equation, in the limit 
in which p = (i.e., p <C mc), can be obtained in the following way: 
In this limit ^ is equivalent to 

dttp = - V ^-(Ty^* , (5) 

where a y is the Pauli matrix. Notice that the complex conjugate of a y is — cr y , so this 
entails 

d t ip* = —j- a y^ > ( 6 ) 



and as a consequence 



,2\2 



W = -[^-)l>. (7) 
The solution must therefore be of the form 

■0m(*) = cos uti/j M (0) + sin ut£ , (8) 

with uj = mc 2 /h, and £ is fixed by consistency to — oy^^O), thus arriving at 

ipM(t) = cos Lut^ M (0) — smuta y ijj* M (0) . (9) 

This is to be compared to the corresponding rest mass frame Dirac equation, which is 
given in this case by 

dti>D = -iuo- z ip D , (10) 
whence the solution directly reads 

V'-d(^) = ex P (—ioJta z ) iPd(0) = cosuttp D (0) — % sin cut cr z ipEi{Q) . (11) 
We can now compare both solutions, giving 

<Mf) = \ IMt) + M-t)\ - \°* Wd{*) - M-t)) ■ (12) 

Accordingly, what we obtain is that the rest dynamics of the Majorana equation 
is quite more involved than the one for Dirac equation. For Dirac, the positive and 
negative energy branches do not mix in this limit, and the only dynamics is a global 
phase associated to the rest energy of the Dirac particle (positive or negative depending 
on the branch). For Majorana, the dynamics contains one term that is like the Dirac 
dynamics, but other three novel terms: in the first one, the temporal dynamics is like in 
the former one, but reversed (evolution backwards in time). In the second one, there is 
a complex conjugation of the spinor, and a spinor flip: in the Dirac language, positive 
and negative energy components would get mixed by this dynamics; in the Majorana 
equation this is more subtle, given that there are not energy eigenstates associated 
to this equation and positive and negative energy branches do not have here a clear 
meaning. In the third one, there is a simultaneous time-reversed dynamics, spin flip, 
and complex conjugation. 
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Figure 1. (%p(t)\a z \tp(t)) for states |^(0)) = (1,0) T (Dirac: solid; Majorana: dashed), 
and \ip(0)) = (l,i) T /^/2 (Dirac: dotted; Majorana: dotted-dashed) . 



Summarizing, what we obtain is that even such a naively simple limit as the rest 
limit contains intriguing features for the case of the Majorana equation: interference 
between forward- and backward-time evolving solutions, between real and complex 
spinor components inside the same Dirac-like positive or negative branch, and between 
spinors associated to these branches. 

In Fig. [l]we plot the average value of a z , (i/;(t)\a z \il;(t)) , for states ^(O)) = (1,0) T 
(Dirac: solid; Majorana: dashed), and \ip(0)) = (l,i) T /\/2 (Dirac: dotted; Majorana: 
dotted-dashed). It can be appreciated that, despite the fact that the Dirac dynamics is 
trivial (the average value does not change with time in either case), the Majorana case 
is not: there are oscillations in (ip(t)\a z \ip(t)) even for a particle at rest, and even in the 
case in which initially only the upper component of the spinor is present. 



This is clearly obtained with Eqs. (|9j) and (11) 

(vz) D (t) = ^ D {t)(j z ^ D {t) = ^ D {0)a z ^ D {0) , (13) 
since a z is proportional to the Dirac Hamiltonian; while for Majorana we have 

A{t)<r g ^ M {t) = cos(2wt)^i / (0) ( x 2 ^ M (0) - sin(2u;t)Im [4/(0)^1/(0)] , (14) 
reproducing the behaviour in Fig. [T} 



3. Quantum simulation in trapped ions 

In order to perform a quantum simulation of the rest limit of the Majorana equation in 
trapped ions, one should point out that the operation ip — y Kip = ip* is unphysical. In 
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order to implement the dynamics of Eq. ([9]), we follow Ref. [TTj : We perform a mapping 
from the two-component Majorana complex spinor onto a four component real spinor, 



Accordingly, the rest-limit Majorana equation in the doubled space takes the form 



This equation, as the one with p ^ 0, is straightforwardly simulated with a four- level 
trapped ion or two two-level ions, as shown in Ref. [T7]. In the rest case considered 
here, no coupling to the motion is needed, and just a two-spin Hamiltonian should 
be necessary In order to measure an observable A in the Majorana space, one 
should measure an observable M* AM in the doubled space, where M = (I il) and 
Mipi = ipM- In order to reproduce the dynamics from Fig. [TJ one should measure 
the four-level population by resonance fluorescence with a cyclic transition, following 
standard electron-shelving techniques [20J. 

4. Conclusions 

We have analyzed the strict nonrelativistic limit of the Majorana equation in which the 
particle is at rest. We have shown that intriguing phenomena appear, even for this rest 
situation, as compared to the Dirac equation: Dirac-like positive energy components 
that turn into negative ones, or real and complex spinor components that interfere in 
a nontrivial way. These novel features can be straightforwardly simulated with two 
trapped ions. 
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Appendix 

From the comparison of the Majorana and Dirac equations it is intuitive that the 
ultrarelativistic limit, that is, m — > 0, will coincide for both. The detailed computation, 
however, is instructive. Both equations presented here are exactly solvable; the 
traditional technique for the Dirac case is to notice that the solutions of the Dirac 
equation are also solutions of the Klein-Gordon equation, and use separation of variables. 




( Vl \ 



(15) 



ihd t ?p4 = —mc 2 a x <8> ayijj^. 



(16) 
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We shall present a rather subtle proof that a similar result holds for the Majorana 
equation; this will allow us to assert an ansatz for the solution of the Majorana equation 
that will provide us with the exact solution. We will then take the limit m — > 0, and 
will derive in which way Majorana solutions tend to Dirac solutions. 
Let us write the Majorana equation as 

ihdt^p = —ihca x d x ip — imc 2 a y ip* , (17) 

and use Fourier transform, 

1>(x,t) = £J -|L=^(t)e^ . (18) 

It follows that 

ihd t ip p = cpa x ip p — im(?o y , ^)*_ v . (19) 

Notice that the complex term couples the p and — p components. 

Directly from the previous expression, by taking the complex conjugate and 
interchanging at the same time p and —p, we infer 

- ihd t ip*_ p = -cpa x ip*_ p - imc 2 o- y ip p . (20) 

In this derivation it is relevant to notice that a y is purely imaginary, so ia y is real. On 
the other hand, we can solve i/j* from the expression above, as 

ip- p = ^o-ydtipp - —o- z i\) p , (21) 

p mc 1 mc 

and on differentiating, we obtain 



-,2 



- ihdti>% = —o- z d t ipp + ^a y d 2 ip p . (22) 

Putting all together, we conclude that the components t/j p obey the Klein-Gordon 
equation 

h 2 d\i) p + (p 2 c 2 + m 2 c 4 ) V P = . (23) 
It follows that ip p (t) will have the solution 

i> p (t) = cos ( Up t) ^p(O) + ^(0) , (24) 

u p 



where u p = y/p 2 c 2 + m 2 c 4 /h, and using the equation above, 

icp mc 



1^(0) = -^vMO) - ^^C P (0) , (25) 



so 



i/j p (t) = cos (u p t) ^ p (0) - Sm ^ pt) [icpa x ^ p (0) + mc 2 a y TP*- p (0)} . (26) 

The limit p ^> mc now makes sense, and gives (together with a condition over the time, 
t < 2h\p\/m 2 c 3 ) 

ip p (t) « exp(-icpta x /k)il) p (0) , (27) 
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which coincides with the aproximate solution in the same limit for the Dirac equation. 
The exact solution for the Dirac equation in all the range of values of p is actually 

tff(O) , (28) 

which allows for direct comparison with the Majorana solution. For example, clearly 
when m = both solutions coincide. 
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